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ABSTRACT 


Three techniques for finding the eigenvalues and eigenfunctions 
are investigated, A typical problem invoives a linear homogeneous 


differential equation with variable coefficients of the form 
P(x) y"(x) + P'(x) y'(x) + wo M(x) 20. (1) 


The functions P(x) and M(x) are functions which are positive, or have 
at most isolated zeroes on the fundamental interval (0,L); w is a 
parameter. Appropriate end conditions are specified so that the 
problem is self-adjoint. 

The three methods are: Newton's method, Stodola's method, and 


the Rayleigh-Ritz method. The metheds cre derived and a computer 


solution by each method is included in the paper. A second problem 
involving Besselts equation of order zero is soived using each method 
and a comparison of the eigenvalues and eigenfunctions is made witn 
tabulated values. 

The results indicate that Newton's method is to be prefezred 


Meu! ly, 








Page 

We SEERITT ER Co) (90@ UEC) gam SRS an Sarai a ks i JE ss igh ee 5 
TI. NEWION'S METHOD - ------------+---+----- 10 

Ae SSSI RUE ON (Ox? INSCINSISee (GNP IEUIRO| SS 10 

B. COMPUTATIONAL ROUTINE -~----.=----~--=--== 13 
TTI, STODOLA'S METHOD - ------------s 4+ - wae 1 
TV.  RAYLBIGH-RITZ METHOD ----.~-----+--------- 22 
V. APPLICATION OF RAYLEIGH-RITZ, STODOLA, AND 

NEWION'S METHODS ----+--++-4++--- a = = = 26 
SEC OCT UCHON Cotten a es oe ee 30 
COMPUTER OUTPUT 1 ----+------+---~-----+----+.-+- 32 
PIGS. (GUIS Sup 2) eS 5 ee 33 
COMPUTER OUIPUT 3 --=-={-+-----+-----4-+--+-+--+---- 3h 
COMPUTER OUIPUT 4 =-----+--+--+-+--+--+-+-+-------= 35 
COMPUTER OUTPUT 5 ------------++-+---+- - 36 
COMPUTER OUTPUT 6 ----~-----+--+-+-+-+-+-+---- 37 
COMPUTER OUTPUT 7 ------------------- -- 38 
COMPUTER PROGRAM 1----~------------------ 39 
COMPUTER PROGRAM 2--------+------------- hy 
COMPUTER PROGRAM 3~----------+------------ 43 
COMPUTER PROGRAM 4 -------------------- - As 
COMPUTER PROGRAM 5--------------------- 46 
COMPUTER PROGRAM 6--------+-----------+---- 48 
COMPUTER PROGRAM 7--------------- eee eee 50 


TABLE OF CONTENTS 








BIBLIOGRAPHY 


INITIAL DISTRIBUTION LIST -------+--+-----+------ 


FORM DD 14773 





IT. INTRODUCTION 


Consicer the linear homogeneous differential equation with 


variable coefficients of the form 


P(x) y (x) + Pt(x) y"(x) + 0° M(x) y(x) = 0, xe(0,1), (1) 


where P(x) is of class gle M(x) is continuous, and both functions are 
positive, or have at most isolated zeroes on the fundamental interval 


(0,L). Alse given are end conditions of the form 


l 
2 


1 
aa y(x, ) + by As (x,) iz 


(2) 


| 
© 


a, y(x,) + by y' (x5) = 


The function y and the parameter w are to be determined. 
An equation of the form (1) may arise as a result of solving a 


partial differential equation of the form 
ne ~ 
n(x) BED - S (p(.) ALbatd) (3) 
ares Ox Ox 
6 t 
by separation of variables with boundary conditions corresponding to 
(2). The solution of (43) is assumed to be of the form 
U(x,t) = X(t) T(*) (4) 


and standard techniques for separation of variables [8] are employed. 
Equation (1) may also occur as the Evler equation [1] of an isopermetric 


problem of the form 





iF : D 
bi P(x) y'(x)* dx = MIN (5) 


subject to the constraint 


= 2 
J M(x) y(x)" dx =1. (6) 


This is equivalent to solving the problem 


L 
J P(x) y' (x) ‘as/f y(x)ax = MIN. (7) 


In this case, the Euler equation is 


ax Bys) - FB, = Ge (PG) yt) ] + ox) ¥(x) = (8) 
where | 
a ae (9) 


f° = P(x) yt (x)* 


(10) 


cel M(x) y(x)*. 


Calculus of variations treats the minimization or maximization 
of functionals such as integrals. This paper considers three techniques 
for obtaining the eigenvalues and eigenfunctions of Eq. (1). These 
techniques are Newton's method, Stodola's method, and the Rayleigh- 
Ritz method (3, 2, 5]. 

Newton's method may be applied to the solution of either the 
isopermetric problem or the equation resulting from separation of 
variables, The principal problem is to obtain a value of the solution 
to the differential equation which will satisfy the end conditions 


imposed by Eq. (2). 








Consider y as a function y(x,”) on those extremals [1] from 


(x, , yy) Let the derivative of y with respect to w be 


, ay (x0) 
S(x,w) a aw ’ (11) 


A correctional equation 


y (x, 9%) ~ y(x,) 


“New ~ “ola ~ S(x, 51) (12) 


is derived in Chapter II and the sequence of steps to obtain the 
solution is discussed there. 


Now consider Eq. (1) written in the self-adjoint form 


- & (PE) yt) | = 0? He) ¥) (13) 


: ; a Ls - 5 = ‘5 Paz 
Subject to the constraints (2). In Stodola's method, consider (13) to 


have the form 
: 2 

Lyswy . (14) 
where the operator L is defined 

Ly(x) = aay me (P@) v'@) (15) 

else dx : 
For the eigenvalues ws and associated eigenfunction Xi, 
; eee 


Taking into account the end conditions, define an inverse operator, 


in then 


it es (17) 
W 





Since differential equation (13) and the boundary conditions (2) 
are self-adjoint, there is an infinite set of orthogonal eigenfunctions, 
Xi, and an arbitrary function, Y,, satisfying the end conditions (2), 


can be expanded in a seriés of eigenfunctions 


Repeated application of the inverse operator C7) tO Yo emphasizes 
the coefficient of X,, and decreases the relative size of the other 
coefficients, The resulting function Y is the corresponding eigen 
function. Further details of Stodola's method are discussed in 
Gaapter Jit. 

The Rayleigh-Ritz method is a procedure for obtaining approximate 
solutions of variational problems. The procedure consists of assuming 
that the desired extremal, y(x) can be approximatea by a linear 


Scmpination of n suitably chosen functions 
y (x) ~ Co, $, (x) + Cy $,,(x) + see + Ce es. : (19) 


The C. are to be determined to effect the desired minimum, Usually 
the 4 i(x) are chosen to meet the boundary conditions for any choice 
of C5. An approximation of y (x) may also be made by spline functions 


of the form 


y (x) SC, +C, x + Cz x + Cy | (x _ 2, ) | + see + Crs | (x - a,)° | 


where the C. are to te determined and the a, are in (0,L). ‘The 


Zi 
function | (2 - %,) | = (x - ie if x >, and is zero for x <a. 





The C's are not independent but must be chosen to satisfy the end 
conditions, In both the polynomial and spline function approximations, 
determination of the C. is accomplished by solving an algebraic 
eigenvalue problem [4]. 

In the following chapters, Newton's, Stodolats and the Rayleigh— 
Ritz mentod are discussed in more detail end computer solutions of 


(1) by the three methods are presented. 





IT. NEWTON'S METHOD 


The problem of minimizing the numerator of Eq. (1.7) with the 
denominator constrained to equal one is again being considered. The 


problem may be expressed in the form 
L 5 L 0 
eS | Gs) ile) sBee l aP cbs = nm (1) 
0 O 
and 


2 a 
en a i ama chee 4 (2) 


x 
Asswne the desired solution y(x) = y (x) has been found. Then 


(3) 


ey 
NO 
ed 
Fla 
Ps 
NO 
ae 
ot 
cS 
NO 
* = 
reo 
Ma 
NO 
j 
>) 


and 
L 1 : ; 
Peas) eke a" 2 J Fy ae a 8 (4) 
O 
A curve satisfying conditions (3) and (4) is called admissable. 


A. DERIVATION OF NECESSARY CONDITION 
The first necessary condition is derived for the case where 
% x ' 
b, =0=b,. Consider replacing y (x) by y (x) + eT] (x), where 


N(x) is an arbitrary function vanishing at 0 and L and having 


piecewise continuous first derivatives. It is necessary that 


L x 
ae J yao ft, ax) WG) ax = 0 (5) 





Or 


: L 

call ! 

ge = J M(x) T (x) dx = 0 
mer all 1 such) thas 

L x 

dJ 1 1 ' 

= 2 J Coe J f', dx) 1 (x) dx = 0 
or 

dJ : : 

ge = J MG T) (x) dx = 0 


fmoeig. (6) and (8), M, (x) is defined as 


e x a 
ul iE 
hae ie 18 eh * 
i) = Gy i” y My ay 5 yt ayer, 
het 


L 
g, (x) = M, (x) - a = M, (x) ~ J M, (x) dsc) ines 


(6) 


(7) 


(8) 


(9) 


(10) 


The above necessary condition may be reexpressed. It is necessary 


ate, G 

ee 

a> =] &(x) 1 (x) ax = 0 

O 

mer ail i} such that 

aI 

ae 4, Pee) a = 0 
since 


L L 
J Deere 1 (x) dx = M. J, N' (x) 0 My 


lav av 0 


Now choose 


M(x) = g(x) + ¢ 6, (x) 


IESE 


(11) 


(12) 


(13) 


(14) 





t 
and choose ¢c so that 7| (x) is admissable: 


L 
dJ , A a 


ie g, (x) #0, Eq. (16) can be solved for c giving 
L L 5 
0 


Now consider 


L 
ss. a J 9, (g(x) +c e, (x) ) dx (17) 
then 
ae L 
= +c =. - J (go) ere, (x) )( go (x) + ¢ £, (x) ) dx (18) 


NO 


f (ql) + ce) azo. 


@) 
Then 
dit 
ne 20 (19) 
unless 
8, (x) + Cc g, (x) =O, (20) 


Hence this is a necessary condition, known as the First Necessary 
Condition, or the integrated form of the Euler equation. Writing 


Eq. (20) as 


AO 1 : | 
fa + C are ~ J ve + Cc f.) dx = (24, (x) + C (=) | (21) 


it is seen that f° and r! enter only on the combination 


ame a (22) 


WA: 





The integral form of Euler's equation [1] 


ak 
mh e = 
Fe a Fy dx = M = CONSTANT (23) 
is obtained from (20) by setting 


Now let 


C= = (25) 


then the Huler Equation 


uw t 
Fytyt y (x) + Foty y (x) + Hee y(x) = i ate, (26) 
becomes 
P(x) y"(x) + PM(x) y"(x) = = wo Mx) y(x) . (27) 


It may be noted that this differential equation is self-adjoint. 


B, COMPUTATIONAL ROUTINE 
It will generally be necessary to generate the curve by forward 


numerical integration. There are two parameters 


Ce cate wr (28) 
and 


>= (0) (29) 


QO 
It 


which are needed to determine a curve which is an approximation to 
the solution; an estimate of yp(0) is needed in order to start 


integration. 





In order to apply Newton's method it is necessary to evaluate 


derivatives with respect to these parameters. Let 


o)peai (ec : 
TL a 3 CG (30) 
1 
and 
0 y(x 
Va 3 C e id 
2 
Then differentiating Eq. (27) with respect to c, gives 
ft t 
P(x) u +P (x)u - C, M(x) u = M(x) y(x) (32) 
Denote the left side of (32) as Lu, defining the second-order 
linear differential operator L. Then the equations for u may be 
written 
Lu =M{x) y({x) oo), 
and 
u(O) = 0 
(34) 
t 
u (0) = 0 
The corresponding equations for v are 
tf | t 
P(x) v +P (x) V - Cc, M(x) v = 0 (35) 
and 
v(0) = Q 
(36) 
{ 
7 (0) = 4 


14 





Oue 


Lve=Q0Q 
0 (37) 
v(0)=1. 


The correctional equation for J is obtained from the approximate 


relations 





el al OJ 
OS ol es, (38) 
=A, - J; 


by J is meant the computed value. The correctional equation for y 


is obtained from the linear approximation 


A y(Z) = ui) Oc, + v(t) de, (39) 


- y(L) 


Il 


where y(L) is the computed value. 

However this problem has a peculiarity that makes it sone 
simpler. Because both integrals are quadratic and homogeneous in the 
pair y(x) and y (x), the principal problem is to find c, so that 


y(L) = Q, The constant Cr is then obtained by substituting the 


Seruyi0on y into J: 
Cy = dd (40) 


Consider then y as a function of a single parameter 


y = y(x, c,) (41) 





The correctional equation reduces to 


Ay = u(L) do, (42) 
= y(L)/x(L) 
Mies = ya) (43) 


In the computational routine, an initial estimate of we is made, 
Then y(x) ard u are obtained by integrating Eq. (33). The correction 
ong we is then obtained from Eq. (43). If a good guess for we is 
made, this method converges rapidly. If a normalized value of y(x) is 
desired then a new solution may be obtained from 


4 


Ly ae 
0 


v(x), = ¥(@)/(J uC) ye)? ax)” (44) 


new 


To tell if the lowest eigenfunction has been obtained, a check 
is made to see if the eigenfunction obtained has a zero in (0,L). 
If it has no zero in (0,L), the. lowest eigenfunction has been 


obtained. 


16 





ITI. STODOLA'S METHOD 


Consider Eq. (1.1) 


& (P(x) y'(x)) = -0® M(x) y(x), (1) 


Subject to the constraints 


Il 
oe) 


a, y(0) +b, y (0) 


(2) 


I 
© 


a, y(L) + by y (1) 


stodola'ts method of finding the lowest eigenvalue and eigenfunction 


will now be developed. The system of Eqs. (1) and (2) is self-adjoint. 


Gh this case there arc an infinite number of eigenvalues Ai = w,* 5 
i 
i= 1, 2, ... and associated eigenfunctions Xo. 
Consider Eq. (1) as having the form 
L y(x) = 0 y(x) (3) 
(oe) ea ie ke 


in which the linear differential operator L is defined by 


b y(x) = h- & (P(@) ¥'()) (4) 


and consider any particular eigenfunction Xi. Then 


ham one o8 
1 
2. A 
L Xie ws Xi (5) 
i eae om Xi 


17 





The inverse operator, ee is now defined, taking into account 


the end conditions (2). Consider 


& (P(x) xi) ue ws M(x) Xi (6) 


and integrate (6) to obtain 
2( . 
P(x) X, = - w; a M(x) Xi dx - c) : ee 


the constant c is to be determined later, The eigenfunction, X., 


may be obtained by dividing Eq. ) by P(x) and integrating again 


= M(x) X; dx 


aout is. af sty dx - a | (8) 


The constants of integration, c and d, mav be obtained from the 


UC... & ol iam inl 
FOENVIVC YT Wee’ a xe 
Col KLIS ale i IS a 


a NE al PN ana (A) Ans 
Ge \-/ ; cet Tél rar] Gane ( !o roma Ata LU e op VE rang these 


equations gives 
L I 
ae M(x) X, dx Ee [ MG) X, ax 


O 
af te, ST - 


a,b - 
Za. | 
la, f 63 p a "2, P(O) 
and 


d=-b, d/a, P(O) (10) 


if a, A705) Ef a, = 0 or P(O) = 0 then c = O and a solution for d 
is obtained. 

The set of operations on the right side in Eqs. (7) and (8) define 
the inverse operator, nl with c and d given by Eqs. (9) and (10). 


These operations give 


18 





2 
eS se L © Gee 


or 


~| 2 
LX; = X, /w, Ge) 


Stodola's method makes use of the fact that an arbitrary piecewise~ 
continuous function which satisfies the end conditions (2) can be 
expanded in a series on the eigenfunctions. The first estimate Yo 
of the first eigenfunction may be chosen as follows. let Yo be the 
first estimate of X.3 choose it so that it has no zeroes inside the 


interval (0,L). Then Yo may be considered to be the following form 


Leer el 


eerie, “yaa nae (13) 


ee 


where a, is not equal to zero. For convenience Yo is normalizeds; the 


4 


lmeciorm of Yo, io lie is defined as 


MAX |Yo(x)| = [|¥o|] . (14) 


Ors Ae i 


and yo | | is Set equal LO ones Now consider 


4 a 
and 
a a 
—m 1 a 
W Ww 
1 2 
2 w. 2m 
apa | er ee) ) 
~ Om oa A, 23 a) C ) 9 ae, ( ) s + eee 
sae ? 
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Fron Eq. (16) it can be seen that the relative size of the components 


DAN Ne 
Koy io eos 18 decreasing by a factor (— - () » eo respectively. 
We Ww 


After a few applications of the operator, glo GO BO sual) eerie 
leading term will) dominate. It is convenient to normelize after 


each iteration; let 
arent Y (17) 


f 
where Ym is the m th approximation to xX, and 


Wm 4 / | 2a (18) 


When this is compared with Eq. (12), where Xx, on the left hand side 


corresponds to ie and X. on the right hand side corresponds to caw 


1 


* * 2 e e 
it is seen that oF 1S approximately 


we = 1/{[Z1| - (19) 


The iteration is stopped when 
[ ly) ~ ees) tee (20) 


where € is some specified small number. The resulting function Y_ 


1S an approximation to Kae 


In the computational routine an initial approximation, Y, to 


X, is made which satisfies the constraints (2) and = is normalized. 


The next approximation to xX, is obtained by forward integration of Kas. 
(7) and (8). The contrants of integration, c and d are determined 
by Eqs. (9) and (10). The value of w 


The new approximation, A is normalized, Eq. (18), tr procedure is 


is approximately /|IZ 1 I- 


20 





repeated. When the norm of the difference between successive iterates, 
a4 Y? is less than some prescribed value the routine is stopped. 
The yalue of w° from Eq. (19) is taken as the lowest eigenvalue and 

Th is the corresponding eigenfunction, The numerical solution of 


Besselts equation of order zero gave a gocd approximation to the 


eigenvalue but a smaller value of the eigenfunction. 


al 





IV. RAYLEIGH-RITZ METHOD 


The Rayleigh-Ritz method has long been used to obtain an approximate 
solution to sigenvalue and eigenfunction problems, This is still 
a useful method if a high speed computer is not available. 

Assume that m functions 2. (x) are given which satisfy the end 
conditions. These are often chosen as polynomials or trigonometric 
functions. The solution then is approximated by a linear combination 


of these 


aa - 
J = 2 a4 § . (x) 9 (1) 
i 


This is substituted into 
Ss ne 2 
Wo” ss i) Pieamy, (35) sax / if M(x) y(x)~ dx = MIN (2) 
@) 0 


and the c's are chosen to minimize (2). This is equivalent to the 


minimization of a quadratic form 


where 
L 
a.. = r Bese) 20.) 6. (eae (3) 


subject to a constraint that 


m m 


3 SM Ge oat (4) 
| ee. 


ee 





where 
L - 
Ds = J M(x) §. (x) (x) obit (5) 


Equation (2) may then be written 





wo EE | MON oO 
CBG 
where 
A= tig 9 La 1s2yeeey my J = 1525-0 om, (7) 
B= Dey 2 La 1y2yeee, my J = 1,2ye00ym, (8) 
and 


m 
ae 


(9) 


p= \ 
CG = ce @ Caos 2¢9e@ Cd 
A = ait 


The problem may be reduced as follows. Let the eigenvalues of 


B be Me > 1 = 1,2,---m and the associated normalized eigenvectors 


ae. Then 
a 


i (10) 
wnere 
T = (ups Uy, eee, U) DIAG (1/dyy 1/dgy oes T/A) (11) 
and the transformation 
where 
= i? 
D= (d,, dng eee, dm) (qs) 


= 





gives 


at (14) 


The desired minimum value of 


cxt=aD PT ATD =100 (15) 


subject to (14) is the minimum eigenvalue of 


E-@T 


+3] 


(16) 
and the corresponding vector D is the associated normalized eigen- 
vector of E. Since the elements of T and D are known, C may be 
evaluated from (4.12). 


The minimization technicgue using the spline function approximation 


y(x) © C, +0, X40 x +c, Rese) Soo G, 


yX+ 05 4,)°) 


(17) 


x43 | 


1s Similar to the polynomial approximation method. In this case 
however, the functions @. (x) do not satisfy the end conditions. 

Hence the constants above are not independent but one cf them is 
dependent on the others for the end conditions to be satisfied. It 
is necessary to have an initial subprogram to eliminate this constant. 


Then the problem is reduced to the same form as pefore 
m-1 m—-1 
Wea coe gene EL ul (18) 
1 og 


i al 
i=] j=1 J 


el 





and 


m—-1 m1 
> t bea Cc. ea (19) 
i=1 j=1 rs 


and solved as before. 


2? 





V. APPLICATION OF RAYLBIGH-RITZ, 


SLODOLA AND NEWTON METHODS 


in Chapters II, III and IV, three methods were developed for 


finding the lowest eigenvalue and eigenfunction for Eq. (1.1). These 


methods were used to obtain numerical solutions of the following 


two problems: Bessel's equation of order zero 


Xy (x) +y (x) + 0° xy(x) =0, xe(0,2-405), 


subject to the constraints 
t 

y(0) = 1, y (2.405) =O, ¥ (0) = 0; 

and an equation of the form 
ff q 2 

(14x) y (x) + y (x) + Ww y(x) 20, = xe(0,1) , 

subject to the constraints 
y(0) = y(1) = 0. 


In this chapter the computation and numerical results 


(1) 


(2) 


(3) 


(4) 


are discussed. 


In the Rayleigh-Ritz polynomial approximation, an approximation 


for the extremal, y(x), satisfying Eq. (3) was made by choosing 


suitable functions 8, (x) and Ayes) to satisfy the constraints in 


Eq. (4). The approximation for y(x) is 
y(x) = C, 8. (x) + Cy (x) 


Cc 


4 x (1-x) + Co a (1-x) 


(5) 





me values of Cas Cos and the lowest eigenvalues were obtained by 
the techniques described in Chapter II. Computer Output 1. shows 
the values of Cy and Cos the lowest eigenvalue and the values of 
y(x) at tenths of an interval. The numerical solution of Eq. (3) 


using the spline function approximation 
2 a ne 
y (x) so, K+ Cy x + Cz | (x=1/3) | + Cy | (x-2/3) | (6) 


is shown in Computer Output 2. and Computer Program 2. This solution 
followed a format similar to the polynomial approximation solution. 
However, a transformation was required to eliminate the dependent 


coefficient. In this case the coefficient c, was taken as dependent 


4 


49 Co and oe The required matrix products, transpositions and 


eigenvalue solutions were evaluated by using subroutine GMPRD, GMIRA, 


ei.C 


and BIGEN respectively, [4]. 
A numerical solution of Eq. (3) by Stodola's method using as a 


first estimate 
y(x) = 4x (1-x) (7) 


to satisfy the constraints in Eq. (4) yielded the values in Computer 
Output 3. The solution of Eq. (3) by Newton's method with an initial 


estimate of 

w’ = 2.0 (8) 
is shown in Computer Output 4. In these two methods the fundamental 
interval (0,1) was sub-divided into one hundred equal sub-intervals 


to perform the second order Runga-Kutta integration routine. if the 


error between the computed value and the boundary value at X equal 


eal 





one was less than one ten thousandth, the program was terminated. 
The programs for Stodola's and Newton's method are shown in Computer 
Program 3. and Computer Program 4. 

The lowest eigenvalues obtained by the above methods are reasonably 
close in value. The maximum difference in values was between Newton's 


method 


w = 3.78634 (9) 


ang. Rayleign-Ritz polynomial approximation, 


Ww = 3.79813 . (10) 


Newton's method converged in three iterations and Stodola's method 


required seven iterations. 


13 


The Rayleigh-Rite, Stodola and Newton methods were used to find a 
numerical solution to Bessel's equation of order zero, Eq. (1). A 
tabulated solution, L7], of Eq. (1) gives the smallest eigenvalue 


as one. The result obtained by the Rayleigh-Ritz metnod was 


= 1.01435 Gap) 


when the function was approximated by quadratic splines 





me). = ©: boy X 40, (x - 0.8)°] +c, (x - 1.6)* | (12) 


The smallest eigenvalue of Eq. (1) obtained by using Stodola's 


method with an initial estimate of 


y(x) = ((1-x)/2.405) (13) 
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w = 1.00000 (14) 


In Newton's method an initial estimate of 


w= 1,4 (15) 


yielded a value 


1.0000 (16) 


= 
i 


for the lowest eigenvalue. 

A comparison of the eigenfunction, y(x), obtained by the three 
methods was made with tabulated values for y(x) of Eq. (1). The 
values obtained by Newton's method et tenths of an interval were 
within five ten-thousancths, by Rayleigh-Ritz were within five 
thousandths and by Stodola's method were within two one-hundredths. 
The best approximation of the lowest eigenfunction was obtained by 
Newton's method. 

The solution of w* obtained by the Rayleigh-Ritz meshod is 
greater than the theoretical vaiue as given in Reference 5. The 
numerical solution of Bessel's equation of order zero supported this 
conclusion. The value obtained for we? was approximately one-thousandth 
larger than the theoretical value. 

The results ot the three methods are shown in Computer Output 
5-, 6., and 7. and the programs are shown in Computer Programs 5., 


ey, and 7. 


2g 





VI. CONCLUSION 


The three methods each yield satisfactory values for the eigen- 
values and the eigenvectors. Generally, Newton's method seems to be 
most satisfactory. It is for the most part straightforward to program. 
The only significant errors that are not apparent are there due to the 
integration routine. It is very flexible, being applicable to many 
different types of problems. It converges rapidly if a good initial 
estimate is made, and computing times seem to be generally small. 
Proof of the exact convergence of Newton's method may be found in 
Reference 4, 

Stodolats method also is quite satisfactory. There seems to be 
x convergence... The imiwial euessaeaa be 
and it will still converge. More iterations than for Newton's 
method seems to be required to get the same accuracy, and there is no 
Simple way to get an estimate of the error. 

The Rayleigh-Ritz method seems to have little to recommené it 
if a good computer is available. It is more difficult to program. 
it has one advantage that no iteration is required. Use of spline 
functions increases the tedium of programming so that they are not 
worthwhile. 

The methods of Stodola and Rayleigh-Ritz were used before the 
advent of large scale computers, though application was rather limited 


compared to problems that can be treated now. 
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in case larger eigenvalues and associated eigenfunctions are 
needed, Newton's method yields these directly if the initial estimate 
of w is changed, The number of zeros of the eigenfunction obtained 
show which eigenvalue and eigenfunction has been obtained. The 
eigenfunctions are automatically orthogonel with a weight factor 
M(x) since the system is self-adjoint. 

If Stodola'ts method is used to obtain higher eigenvalues and 
eigenfunctions it 1s necessary to adjoin a further condition that 
each new eigenfunction be orthogonal to all of the preceeding. ‘This 
increases the tedium of programming and decreases the accuracy, 
since preceding values are only approximate. 

in the Rayleigh-Ritz method, an estimate of other eigenvalues 
and eigenfunctions is obtained from the ovsher eigenvectors of the 
matrix E of chapter IV. ‘The eigenvalues are toc high unless the 
eigenfunction can be expressed as a finite sum of the functions 


$. (x). 
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